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Abstract 

Using the factorization approach and considering the contributions of the current- 
current, QCD penguin and electroweak penguin operators at the leading approxima- 
tion, the decay amplitudes and decay widths of Bq — > Dq{D*)P and Bg — > Dq{D*)V 
transitions, where q = u,d,s and P and V are pseudoscalar and vector mesons, are 
calculated in terms of the transition form factors of the Bg — > Dg and Bg — > D* . Hav- 
ing computed those form factors in three-point QCD sum rules, the branching fraction 
for these decays are also evaluated. A comparison of our results with the predictions 
of the perturbative QCD as well as the existing experimental data is presented. 
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1 Introduction 



With the chances that a lot of Bq mesons will be produced in B factories [1,2], it would 
be possible to check the two-body non-leptonic charmed decay modes Bq ^ Dq{D*)P and 
Bq — > Dq{D'^)V . Analyzing of such type decays could give valuable information about the 
origin of the CP violation, hadronic flavor changing neutral currents, test of the standard 
model (SM), constraints on new physic parameters as well as strong interactions among the 
participating particles which provides valuable tests of the QCD factorization framework. 

Theoretically, analyzing of the two-body B-decay amphtudes have been started using the 
framework of so called "naive factorization" [3-7]. This method for some decay channels 
is replaced by QCD factorization [8,9] since it could not predict direct CP asymmetries 
in those decay modes. First, the QCD factorization approach had been applied for the 
simplest charmless S — >• tttt and B — > t^K decays [8, 10-13] then extended to the vector 
and exotic mesons in final states [14-17] and r] or rj' with a pseudoscalar or vector kaon 
[18]. In [19-21], decay modes of Bs meson are discussed. A comprehensive study of the 
exclusive hadronic B-meson decays into the final states containing two pseudoscalar mesons 
(PP) or a pseudoscalar and a vector meson (PV ) is discussed in [22]. The Charmless anti- 
cs — > VV decays has also been analyzed in QCD factorization in [23] . The hard-scattering 
kernels relevant to the negative-helicity decay amplitude in B decays to two vector mesons 
are calculated in [24] in the same framework. The two-body hadronic decays of B mesons 
into pseudoscalar and axial vector mesons have been studied within the framework of QCD 
factorization in [25]. A detailed study of charmless two-body B decays into final states 
involving two vector mesons (V V ) or two axial- vector mesons (AA) or one vector and one 
axial-vector meson (V A) has also been done within the framework of QCD factorization 
in [26]. Considering the contributions of both current-current and penguin operators, the 
anii:)litudes and branching ratios are recently estimated at the leading approximation for 
B^ ^ B*P,BV in [27]. 

In the present work, taking into account the contributions of the current-current, QCD 
penguin and electroweak penguin operators at the leading approximation, we describe the 
charmed decays Bg — > Dq{D*)P and Bq — > Dq{D*)V in the framework of the QCD factor- 
ization method. First, using the factorization method, we calculate the decay amplitudes 
and decay widths of these decays in terms of the transition form factors of the Bq Dq and 
Bq — > D*. Having calculated these transition form factors in the framework of the QCD 
sum rules in our previous works in [28, 29], we calculate the branching ratio of these decays. 
In order to estimate the approximate branching ratios and to have a sense of the order 
of amplitudes, we make a rough approximation, i.e. at the leading order of a^. Within 
this approximation, the hard-scattering kernel functions become very simple and equal to 
unity [27]. In this approximation, the long-distance interactions between the P{V) and 
Bq — Dq{D*) system could be neglected. The higher order ag corrections might not be 
small, but calculation of these contributions is not as easy as the light systems in final state 
and needs much more efforts. Hence, for obtaining the exact results on the considered tran- 
sitions, those effects should be encountered in the future works. There are several methods 
in which such type contributions can be studied, QCD Factorization [8-10, 18,22,24], Per- 
turbative QCD [30] and Soft-CoUinear Effective Theory [31,32]. For more detail analysis 
of NNLO corrections to S — > light — light systems and higher order QCD corrections to 
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the charmless B decays see also [33-38]. Note that, some of the considered decays in this 
paper have been analyzed in the framework of the perturbative QCD in (PQCD) [30] and 
for some of them, we have some experimental data [39]. 

The outline of the paper is as follows: In section 2, we calculate the decay amplitudes 
and decay widths for Bq Dq{D*)P and Bq ^ Dq(D*)V transitions. Finally, section 3 is 
devoted to the numerical analysis, a comparison of our results with the predictions of the 
PQCD as well as the existing experimental data and discussion. 

2 Decay amplitudes and decay widths 

In the present section, we study the decay amphtudes and decay widths for Bq — > Dq{D*)P 
and Bq Dq{D*)V decays, where q = u,d,s , P = tt, K, Dg> and V = K*, D*, {q = d,s). 
At the quark level, the effective Hamiltonian for Bq — > Dq{D*)n{K, K*) is given by 

Heff = ^{V,,V\{C,0^,+C,0^,)}. (1) 

Here O" and O2 are quark operators and are given by 

Oi = {qlui)v-A{cjbj)v-A, = {qluj)v-A{cjbi)v-A, (2) 

where q = d,s and {qiq2)v±A = ^i7'^(l ± 75)?2- However, the effective Hamiltonian for 
Bq ^ Dq{D*) Dg' and Bq ^ Dq{D*) D*, at the quark level can be written as 

Heff = -^^V,i>V:g,{C,0t + C20^,) + Vt,V*,Y.Cn0ny (3) 

Here 0„ are quark operators and are given by 

Oi = {qici)v-A{cjbj)v-A, = {q'iCj)v-A{cjbi)v-A, 

0^(5) = {(libi)v-AJ2q{qjQj)v-{+)A, 04(6) = {qibj)v-AJ2q{qj<li)v-{+)A, (4) 

^7(9) = i(libi)v-AT,qle.q{qjqj)v+{-)A, 08(10) = {qibj)v-AT,qleq{qjqi)v+{-)A, 

where J2q sums over q — u,d, c, s, b and i and j are color indices. The operators Oi 
and O2 are called the current-current operators, O3...OQ are QCD penguin operators and 
O7...O10 are called the electroweak penguin operators. 

The Wilson coefficients C„ have been calculated in different schemes [40-43]. In this 
paper we will use consistently the naive dimensional regularization (NDR) scheme. The 
values of C„ at nif, with the next-to-leading order (NLO) QCD corrections are given 
by [40-43] 
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The decay width and the branching ratio of the nonleptonic process Bq — > Dq(D*q)M, 
where M stands for the P or V mesons, is given by: 



Br{Bq^Dq{D\)M) = TBT{Bq^Dq{D\)M) 

where, A(a;, y, z) = + y"^ + — 2xy — 2xz — 2yz is usual triangle function. 



(6) 



To obtain the decay width, we should calculate the amplitude A. This amplitude is ob- 
tained using the factorization method and the definition of the related matrix elements in 



terms of form factors for the Bq — > Dq and Bq — > D* weak transitions as: 



< Dq{p')\n,{l - ^,)b\Bq(p) (p + pXf:^-''\q') + (p-p')^f^^''''^(q% (7) 
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< Dl{p',e) I 07^756 I Bq{p) > 



mn* + TUB, 



+ 



rriD* + niB, 



(9) 



where is transferred momentum square, q^ — {p — p'Yi and p and p' are momentum of 
the initial and final meson states, respectively. 

We obtain the A as following: 

for Bq DqP (P = TT, K) and Bq DqDq,: 



^^d(s)-*-f^d(s)-P 
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for Bq DqK* and Bq DqD*, : 
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.Bq-*D*D*, .Gf , 

V2 " " 



q ^ q q q ^ ^ 

X [KbV;;, ai - FtftV;;, (04 + a^o)], (16) 



where 



= --J^ , (17) 



In the above expressions, the £*, e , Ek stand for the polarization of the D*, and K* 
mesons, respectively. The quantities Oj, are given in terms of the coefficient Cj, 

1 1 

= Cj + — Ci+i (z = odrf); = + — Ci_i (z = et;en), (18) 

where i runs from i = 1, 10 and 1^^ is number of color in QCD. In the above equation, 
the a\ and are both related to the coefficients C1.2, which are very large comparing 
with the other Wilson coefficients, but we will keep all coefficients to get ride of further 
approximation. 

Now we can calculate the decay widths for — > Dq{D*)P and Bg — > Dq{D*)V decays. 
The explicit expressions for decay widths are given as follow: 
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Table 1: Values of the masses in MeV [39]. 



/.[39] 


hm 


/d±[39] 


/d..[39] 


/dJ39] 


/x*[44] 


130.7 ±0.46 


159.8 ± 1.84 


222.6 ± 19.5 


222.6 ± 19.5 


294 ±27 


217 ±5 


/^.±[39] 




/i.%[45] 


/b±[39] 


/bo [39] 


/bJ46] 


230 ± 20 


230 ± 20 


266 ± 32 


176 ± 42 


176 ± 42 


206 ± 10 



Table 2: Values of the decay constants in MeV. 

3 Numerical analysis 

This section encompasses our numerical analysis, comparison of our results with the predic- 
tions of the PQCD as well as the existing experimental data and discussion. The expressions 
of the amplitudes and decay widths depict that the main input parameters entering the 
expressions are Wilson coefficients presented in the section 2, elements of the CKM matrix, 
leptonic decay constants, Borel parameters and M| as well as the continuum thresholds 
So and Sq [28, 29]. In further numerical analysis, we choose the numerical values as presented 
in the Tables 1, 2 and 3. The Borel mass squares Mf and M| and continuum thresholds 
So and Sq are auxiliary parameters, hence the physical quantities should be independent 
of them. The parameters So and Sq, are determined from the conditions that guarantees 
the sum rules for form factors to have the best stability in the allowed Mf and M| region. 
The working regions for Mf and as well as the values for continuum thresholds are 
determined in [28,29]. Here, we choose the values so = 35 ± 5 GeV^, Sq = 7 ± 1 GeV^, 
Ml — 17.0 ± 2.5 GeV'^, Mf = 7 ± 1 GeV^ from those working values for auxihary parame- 
ters. The values of the form factors f± and /o,i,2,3 at different values of which we need 
in the expressions for decay widths are presented in Tables 4 and 5, respectively. Using 
the expressions for total decay widths, the values of branching fractions for Bg DgP, 
Bq — > D*P, Bq — > DqV and Bq — > D*V are found. We depict the values of the branching 
ratios in Tables 6, 7, 8 and 9. Here, we should stress that, as we mentioned before, our 
results depicted in the Tables are approximate results since we considered the observable 
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Table 3: Values of the elements of the CKM matrix [39]. 
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only at the leading order of as- To obtain more exact results the higher order ag correc- 
tions should be taken into account. However, the presented uncertainties in the results are 
belong to the uncertainties in the values of the input parameters as well as variations in 
form factors which are related to the errors in determination of the auxiliary parameters 
namely, Borell mass parameters and M| and continuum thresholds sq and s'q. These 
Tables also include a comparison of our results with the existing predictions of the PQCD 
as well as the experimental data. From these Tables, we see a good consistency among two 
non-perturbative approaches and the experiment in order of magnitude. In many cases, the 
presented results out of order of magnitude from three approaches coincide especially, when 
we consider the uncertainties in the results. The best consistency between our results and 
predictions of the PQCD is related to the 5° D*^K^ transition and S° D/K*(892)'^ 
transition shows the biggest discrepancy between two methods. Our central value predic- 
tion on — > D*^K^ is approximately the same as the experimental result, however, the 
central experimental result on the branching ratio of — > Dl'^Df^ depicts a big discrep- 
ancy comparing that of our prediction. The presented predictions from PQCD are related 
to the charm- charmless cases in the final states and we have no predictions on the charm- 
charm cases from this approach. In this approach, the wave functions of the participating 
mesons, which are available with higher order corrections, have been used to calculate the 
amphtudes [30] . Therefor, over all agreement between our results and predictions of PQCD 
for charm-light cases in the final state and the experimental data for both charm-light and 
charm-charm cases, could be considered as a good test of the QCD factorization at leading 
order of a^. for related transitions. However, for exact comparison, much more efforts are 
needed in the future works, which may include the higher order corrections. Our results of 
some decay modes which have not been measured in the experiment can be tested in the 
future experiments at LHCb and other B factories. 

In conclusion, using the QCD factorization approach and taking into account the con- 
tributions of the current-current, QCD penguin and the electroweak penguin operators at 
the leading approximation, the decay amplitudes and decay widths of B^ — > Dq{D*)P and 
Bq — > Dq{D*)V transitions were calculated in terms of the transition form factors of the 
Bq — > Dq(D*). Having computed those form factors in the framework of the three-point 
QCD sum rules in our previous works, the branching fraction for these decays were also 
evaluated. A comparison of our results with the predictions of the perturbative QCD as 
well as the existing experimental data was presented. Our results are over all in a good 
agreement with the predictions of the PQCD and the existing experimental data. Our 
predictions on some transitions, which have no experimental data can be checked by fu- 
ture experiments at LHCb or other B factories. To get more exact results from the QCD 
factorization method, higher order as corrections should be considered in the future works. 
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-0.20 ± 0.06 


-0.21 ± 0.05 


-0.22 ± 0.06 


-0.38 ± 0.10 


-0.44 ± 0.11 


-0.49 ± 0.12 


-0.69 ± 0.17 


/f- 




0.58 ±0.16 


0.59 ± 0.15 


0.63 ± 0.17 


0.86 ± 0.22 


0.92 ± 0.22 


0.95 ±0.23 


1.08 ± 0.27 






-0.20 ± 0.06 


-0.21 ± 0.05 


-0.22 ± 0.06 


-0.37 ± 0.10 


-0.44 ± 0.11 


-0.49 ± 0.12 


-0.69 ± 0.17 


/+' (r) 


0.26 ± 0.06 


0.27 ± 0.06 


0.28 ± 0.07 


0.35 ± 0.09 


0.38 ± 0.10 


0.39 ±0.10 


0.42 ± 0.11 






-0.11 ± 0.03 


-0.12 ± 0.03 


-0.13 ± 0.03 


-0.15 ± 0.04 


-0.16 ± 0.04 


-0.17 ± 0.05 


-0.18 ± 0.05 



Table 4: The values of form factors f± at different values of q^. 
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'^ii'*(892)± 






m^,± 




fB+ 

Jo 




0.76 ±0.19 


0.78 ± 0.19 


0.80 ± 0.20 


0.97 ± 0.24 


1.01 ± 0.25 


1.09 ±0.26 


1.13 ± 0.27 


fr 




0.62 ±0.16 


0.63 ± 0.15 


0.67 ± 0.16 


0.98 ± 0.24 


1.08 ± 0.26 


1.14 ± 0.27 


1.27 ±0.29 


fB+ 
J2 




0.90 ± 0.22 


0.96 ± 0.23 


0.99 ± 0.24 


1.50 ± 0.38 


1.60 ± 0.40 


1.82 ± 0.46 


2.00 ± 0.60 


fB+ 
J3 




-1.51 ± 0.38 


-1.62 ± 0.40 


-1.65 ± 0.40 


-2.01 ± 0.50 


-2.30 ± 0.55 


-2.50 ± 0.60 


-2.65 ± 0.61 




0.76 ± 0.19 


0.78 ± 0.19 


0.81 ± 0.20 


0.97 ± 0.24 


1.02 ± 0.25 


1.10 ± 0.26 


1.13 ± 0.27 




-^'\q^) 


0.61 ± 0.15 


0.63 ± 0.15 


0.66 ± 0.16 


0.98 ± 0.24 


1.07 ± 0.26 


1.14 ± 0.27 


1.27 ± 0.29 


fBO- 
J2 


-^'\q^) 


0.90 ±0.22 


0.95 ± 0.23 


0.99 ± 0.24 


1.51 ± 0.38 


1.61 ± 0.40 


1.82 ± 0.46 


2.01 ± 0.60 


fBO- 
J3 




-1.62 ± 0.38 


-1.62 ± 0.40 


-1.66 ± 0.40 


-2.01 ± 0.50 


-2.31 ± 0.55 


-2.51 ± 0.60 


-2.66 ± 0.61 


Jo 




0.38 ± 0.10 


0.40 ± 0.11 


0.41 ± 0.11 


0.58 ± 0.15 


0.62 ± 0.16 


0.67 ± 0.17 


0.70 ± 0.18 




0.33 ± 0.08 


0.36 ± 0.08 


0.40 ± 0.11 


0.67 ± 0.17 


0.71 ± 0.17 


0.74 ± 0.18 


0.76 ± 0.18 




-^'^\q^) 


0.43 ±0.11 


0.47 ± 0.12 


0.50 ± 0.13 


0.81 ± 0.21 


0.85 ± 0.22 


0.88 ± 0.23 


0.91 ± 0.23 




-^'^\q^) 


-0.67 ± 0.17 


-0.69 ± 0.17 


-0.72 ± 0.18 


-1.29 ± 0.32 


-1.42 ± 0.34 


-1.49 ± 0.36 


-1.66 ± 0.36 



Table 5: The values of form factors /o,i,2,3 at different values of q^. 
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DgP 



present work 



PQCD [30] 



Exp 



B^- 
B2- 



B2- 



5.95 ± 1.95) X 10-=^ S.lltiiKilali x IQ-^ 

4.31 ± 1.52) X 10-^ 4.ool?:ilg:ilgi^ X 10-4 

3.44 ± 1.22) X 10-4 
2.03 ± 0.85) X 10-2 

5.69 ± 1.70) X 10-3 2.69ll:I?l°:il°:[ii X 10-3 

3.53 ± 1.23) X 10-4 2.43^foi-a7"i-ao7 x 10-^ 

2.87 ± 0.89) X 10-4 

8.88 ± 2.82) X 10-3 

1.42 ± 0.57) X 10-3 2.13liitlg:KS6 x 10-3 



1.03 ±0.51) X 10-4 
1.20 ±0.73) X 10-4 
2.17 ±0.82) X 10-3 



-1 71+0.92+0.58+0.05 ^ i n-4 
-'-■'-'--0.65-0.55-0.05 ^ -'-'-' 



(4.92 ±0.20) X 10-3 
(4.08 ±0.24) X 10-4 
(4.80 ± 1.00) X 10-4 

(1.09 ±0.27) % 
(3.40 ±0.90) X 10-3 
(2.00 ±0.60) X 10-4 
(1.90 ±0.60) X 10-4 
(6.50 ±2.10) X 10-3 
(3.80 ±0.30) X 10-3 



Table 6: Values for the branching ratio oi Bq ^ DqP. 



14 



B 



d;p 



present work 



PQCD 



Exp [39] 



5° 



DfDs^ 



(4.89 ± 1.52) X 10-^ 
(3.38 ± 1.04) X 10-4 
(2.57 ±0.88) X 10-4 
(11.03 ±2.91) X 10-3 
(3.45 ± 1.75) X 10-3 
(2.08 ±0.68) X 10-4 
(3.14 ± 1.46) X 10-4 
(8.69 ±2.88) X 10-3 
(2.11 ±0.73) X 10-3 
(1.59 ±0.67) X 10-4 
(0.30 ±0.11) X 10-4 
(2.54 ±0.57) X 10-3 



r r,^+2.92+0.44+0.15 ^ 1 r,-3 



q pr,+2.33+0.62+0.12 -, r,-4 
O.UU_]^ g2_o.92-0.12 ^ 



9 «n+l-''3+0.53+0.07 y 1^-3 
Z.OU_i i4_o.70-0.07 ^ -'-U 



9 07+1.52+0.62+0.07 ^ -1 r, 
'^•"J' -0.99-0.69-0.07 ^ 



-4 



9/19+1.12+0.78+0.07 ^ 1 n- 
^•^^-0.72-0.77-0.07 ^ 



1 Rc+0.90+0.56+0.05 ^ -1 n-4 
1.0d_o.63-0.53-0.05 ^ 



(4.60 ± 0.40) X 10-3 
(3.70 ± 0.40) X 10-4 

(10.00 ±4.00) X 10-3 
(2.76 ±0.21) X 10-3 
(2.14 ±0.20) X 10-4 
(9.30 ± 1.50) X 10-4 
(8.80 ± 1.60) X 10-3 



Table 7: Values for the branching ratio oi Bq ^ D*P. 
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B,i D,y presont Avork PQCD [30] Exp [39] 

5± ^ D^K*{m2f (2.90 ± 0.88) x 10"^ 6.49l|il$i^lg:^g x 10"^ (6.30 ± 0.80) x lO'^ 

5± ^ D^D*^ (5.61 ± 1.88) X 10-^ - (4.60 ± 0.90) x 10"^ 

D'^D*/ (7.01 ± 2.09) X 10-3 - (7.20 ± 2.60) x lO'^ 

B° D±K*(892)^ (3.20 ± 1.15) x 10"^ 4.07j:?|^i;:?flg:ii X 10"^ (4.50 ± 0.70) x 10"^ 

S° D^D*^ (8.18 ± 2.84) x 10"^ 

B^ ^ D^D*^ (9.23 ± 2.67) X 10-3 - (8.60 ± 3.40) x 10-^ 

B° ^ Ds'^K*(892f (0.50 ± 0.22) x 10"^ 3.02^i:?^l:g|?ig:ig X 10-4 

B^^^Ds^D*^ (1.07 ± 0.59) X 10-4 

B° ^ Ds^Dl"^ (2.62 ± 0.93) x 10"^ 



Table 8: Values for the branching ratio oi Bq ^ DgV . 
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present work 



PQCD[30] 



Exp[39] 



5± ^ D*°K*{892f 



D* D 



n*± 



(5.07 ± 2.61) X 10-^ 6.82li^^ti:il°Jl x 10"^ 
(0.11 ± 0.07) X 10-2 

(6.85 ± 2.98) X 10-2 



(8.30 ± 1.50) X 10-^ 
< 1.1 % 
2.20 ± 0.70 % 



5° ^ D*^i^:*(890)^ (3.55 ± 1.25) x 10"^ 4.88l|Jilli?lgii x 10"^ (3.30 ± 0.60) x 10 



-4 



50 _^ D*^D*^ 



50 _^ 0*^01"^ 
_^ D;±X*(890)^ 

B^ Dl^DV^ 



(8.78 ± 2.50) X 10-^ 
(8.17 ± 2.93) X 10-2 

^1.63 ± 0.86) X 10-4 3.47l}:ili:°^l°:ll X 10-4 
(6.76 ± 2.69) X 10-4 
(2.77 ± 0.76) X 10-2 



(8.30 ± 1.01) X 10-4 
(1.79 ±0.16) % 



(231?^) % 



Table 9: Values for the branching ratio ol Bq ^ D*V . 
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